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Comment on "Coherence measures for heralded single-photon sources" 

Stefano Bettelli 
Romanogasse 29/15, 1200 Wien, Osterreich 

A recent Brief Report [Phys. Rev. A 79, 035801 (2009)] investigates two figures of merit for her- 
alded photon sources based on spontaneous parametric down- conversion with a continuous pump, 
namely the time-averaged temporal coherence between the signal and idler beams and the time- 
averaged conditioned second-order degree of coherence for photons in the heralded arm. However, 
contrary to what is claimed, no one-arm second-order effect is actually measured. 



PACS numbers: 42.50.Dv, 03.67.Dd, 42.50.Ar, 42.65.Lm 
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The phenomenon of spontaneous parametric down- 
conversion (SPDC) consists in the splitting of a pump 
photon, in a non-linear crystal, into two modes (arms). 
A detection in one arm (idler) implies the presence of a 
photon in the other (signal): the signal photon can thus 
be heralded in a time window of width Td, set by the 
indetermination on the detection time. The probability 
that a second signal photon shows up in the window is, 
for sufficiently large windows, essentially Rt^, where R 
is the pair-generation rate, which can be made as low as 
desired. For this reason SPDC-based sources are often 
used to emulate single-photon sources. 

The apparent Poissonian statistics in one arm of the 
SPDC source [1] observed in the heralded time window 
is determined by the fact that Td is often much larger 
than At, the signal coherence time, typically in the range 
O.lps to lps; only when Td ~ At the thermal nature of 
SPDC radiation shows up, and the probability of a second 
photon in the window approaches 2i?T<j. Experiments on 
SPDC intrinsic one-arm statistics in single-photon regime 
are difficult to perform when sources are operated with a 
continuous pump, since in this case the jitter Td of single- 
photon detectors is not better than ~ 10 2 ps. Thermal 
statistics for one arm of a SPDC source operated in this 
way was demonstrated only recently 

A recent Brief Report by Bocquillon et al. Q (see also 
the expanded version by Razavi et al. [3] ) , concerning the 
characterisation of heralded SPDC-based photon sources 
with a continuous pump, claims ". . . [to carve] into the 
theoretical aspects of [measuring] the temporal correlation 
between the signal and idler beams, and . . . the second- 
order degree of coherence for the heralded signal photons '.' 
More specifically, since the time indetermination of the 
detection apparatus is orders of magnitude larger than 
the coherence times of the source, the authors develop a 
theory connecting the values of the parameters intrinsic 
to the investigated system to the actually observed aver- 
ages, and test the theory experimentally. In particular, in 
the second part of the report, the authors investigate the 
detection of two photons in the signal arm conditioned 
on the detection of one photon in the idler arm; clearly 
this measurement is strictly related to that described in 
[l|, and definitely not less difficult. However, contrary to 
what is claimed, no one-arm second-order effect is actu- 
ally measured in as shown in the following. 



Preliminaries 

The authors of [2| assume a theoretical model for 
a low-gain frequency-degenerate type-II SPDC process 
predicting that all auto- and cross-correlations of the 
source arms can be expressed by means of the two real 
even functions R(t) and C(r) shown in Fig[TJ In prac- 
tice R(t) is the first-order coherence of both arms, and 
R(0) = R ~ 43MHz (i? S PDC in the original text) is 
the pair-generation rate. With the help of the Gaus- 
sian moment-factoring theorem the temporal correlation 
9si (t) of the signal and the idler, Fig|2 is shown to be 
sharply peaked for |r| < At/2 (the parameter At ~ 0.3ps 
plays the role of a signal-idler coherence time): 
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This sharp peak is smeared by detector response times 
and by the finite coincidence window (implemented in 
software), respectively characterised by Td and r c (r CO j n in 
the original text). The authors of [2| model these effects 
by means of convolution functions. For the purposes of 
this Comment it is sufficient to consider r c 3> r^; in this 
case the combined response function of the hardware- 
software setup is equivalent to a moving-window aver- 
age with size 2r c and transition regions as large as 2ra. 
The time-averaged second-order coherence g si (t), pro- 
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FIG. 1: The signal-signal or idler-idler autocorrelation func- 
tion R(t), and the signal-idler cross-correlation function C(r) 
for the SPDC process analysed in || . R is the pair-generation 
rate, while 1/At is the bandwidth. Vertical scales are very 
different since in the actual experiment RAt ~ 1.4 ■ 10 -5 . 
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FIG. 2: The second-order coherence function between signal 
and idler, g^ (t), and its time average, g^ (r). The convo- 
lution function that transforms the former into the latter is 
represented on top. Since X = (2_Rr c ) _1 ~ 10 1 , the smeared 
signal-idler correlation peak, though strongly reduced with 
respect to the ideal one, is still rather simple to measure. 



portional to the measured rate of signal-idler coincidences 
A^j(r), will be equal to one for |r| > t c + r^, and to a 
constant when |r| < r c — Td, i.e., when the large peak of 
( T ) ^ a ^ s m t ne central part of the convolution window, 
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where X — (2Rt c )~ 1 is a signal-idler coincidence factor. 
Since R ~ 43MHz and 2r c < 3ns, the additional term in 
EqJ2]is still of the order of 10 1 , and therefore the signal- 
idler temporal correlation is relatively easily to demon- 
strate. Of course, this signal-idler peak could be made 
arbitrary large by reducing the rate R (in practice, by 
reducing the intensity of the pump, as in 0, Fig.4]). At 
the same time, given that both time scales of the experi- 
mental convolution function, Td ~ 350ps and r c > 0.5ns, 
are several orders of magnitude larger than the coher- 
ence time At ~ 0.3ps, the experiment is only sensitive 
to the integral of the peak, and not to a specific peak 
model. The behaviour of the average g si (r) is illustrated 
in Figl^l and agrees with 0, Fig. 2] (also quantitatively, 
but for the curve at 2r c = 0.78ns for which Td<^ r c ). 

The authors of the Brief Report claim that good agree- 
ment between theory and experimental data proves that 
their model is particularly accurate, but indeed it is only 
the dependence on r c that is truly confirmed, and repeat- 
ing data analysis with different values is just a test that 
the software is working correctly. 



ditioned on idler-photon detection, defined as 
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where the (. . . ) pm averages are performed on post- 
measurement states (i.e., after detection of an idler 
photon). By using the operator identity (Y) pm = 
(E}(ti)YEi(ti)) ■ (Ejit^EiiU)}- 1 , Eqlis reduced to 
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where P SS i is the ideal triple-coincidence rate that, again 
with the help of the Gaussian moment-factoring theorem, 
is proven to be equal to 

P S si(ti,t 2 ,ti) = 2C(t 1 ~t l )C(t 2 -t l )R(t 1 -t 2 ) + 
+ R[R 2 + R 2 (h-t 2 ) + C^h-U) + C 2 (t 2 -U)] . (5) 

Actual measurements are performed with the particular 
case t\ = ti and t 2 = ti + r in mind, so that the relevant 
formula for the numerator of Eq[4] is 

P ssl (r) = P ssl (ti ,U + t, U) = (6) 
= 2C(0)C(t)R(t)+R [R 2 + R 2 (t) + C 2 (Q) + C 2 (t)] . 

This expression cannot be simplified further, but its be- 
haviour is easy to understand. For delays larger than 
the coherence time all non-constant terms die out, and 
P ssi -> R 3 + RC 2 (0) = R 3 g[f(0); for |r| < At the shape 
is peaked, and its maximum value does not exceed four 
times its large-delay value, as shown in FigJ3] The de- 
nominator of Eq[4l i.e. R 3 g^ 2 \o) g^(r), is almost al- 
ways R 3 g^\o), but in the narrow region [— At/2, At/2] 

it grows several orders of magnitude to reach R 3 [g]?' (0)] 2 , 
as can be seen in Figl^J Therefore the function of interest, 
5c (t) = g c (ti, ti + r\ti), represented in FigGH has unit 
value for |r| > At, and is O(RAt) ~ 10~ 5 for \t\ < At/2. 

Experimentally measuring a time-averaged version of 
this well is far from trivial, since the size of the well is 
only O(At), and the typical convolution spread due to 
detector jitter and finite coincidence windows is r c , so 
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Measurement of the conditioned second-order 
coherence 



In the second part of the Brief Report the authors 
analyse the statistics of signal-signal coincidences con- 



FIG. 3: The triple-coincidence rate P S si(t), and the second- 



order correlation function g^ (r) for an idler photon arriving 
at time ti and two signal photons arriving at ti and ti + r, 
as defined by 0, Eq.9]. The temporal correlation function is 
very small (~ 10~ 5 ) in the region r 6 [—At/2, At/2]. 
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that the effect is in the best case 0(At/Td) ~ 10 -3 , and 
even smaiier with the coincidence windows used in the 
article. This is an effect of the same order of magnitude 
as that measured in [l[ on an almost identical physical 
system. Instead, Fig. 3 of [2J shows a large well extending 
over Ins or more. The reason for this strange result is 
that the authors actually plot the "time-averaged second- 

(2) 

order correlation function" cjc (t) defined as 



N sl {Q)-N st { T y 



(7) 



where N S i is the already mentioned signal-idler coinci- 
dence rate, and N SS i is the triple signal-signal-idler coin- 
cidence rate, proportional to a smeared P SS i . It is impor- 

_(2) 

tant to understand that even if g c is studied for t\ = ti, 
the observable N ss { is the outcome of a triple-coincidence 
experiment, and its value is linked to the full three-time 
structure of P SS i, shown in FigJU Minor details apart, 
Pssi{ti,t2,U) is characterised by three levels: the acci- 
dental triple-coincidence plateau at R 3 , the signal-idler 
coincidence ridges at R 3 + RC 2 (0) (when |* x — **| < At/2 
or |ta — ij| < At/2), and a more complicated central peak 
(when both ti, t 2 G [U - At/2, t t + At/2]) where the cor- 
relation grows even more, up to a factor four. 

The experimental rate N ssi is obtained by convoluting 
P ss i with the response functions of three detectors and 
two coincidence windows. The effect can be modelled as 
a moving average with a square bidimensional averaging 
window of area (2r c ) 2 and transition regions of size 2td- 
The value of N SS i will then be significantly different from 
R 3 only when the averaging window is close to one of the 
ridges of FigJU due to translation invariance the result 
is R 3 [l + X] as before. But when the averaging window 
sits at ti = t% = ti both ridges participate to the average 
and their contribution is twice as large: R 3 [1 + 2X}. Note 
that the complicated peak at the ridge junction is almost 
irrelevant since its contribution relative to R 3 X is only 
0(At/r c ) ~ 10~ 3 , due to its height being limited and its 
base being only At x At, and not At x r c as for the part 
of each ridge within the averaging window. 

The conclusion of the previous discussion is that the 
expected shape of the coincidence functions [5(] are 




FIG. 4: The ideal triple-coincidence rate P S si(ti,t2,ti) (only 
most relevant structures shown). The ridges running parallel 
to the axes ti = ti and t-z — ti are given by signal-idler correla- 
tion. When both ti and t% are very close to ti correlations are 
increased by a factor at most four. P sa i is symmetric under 
ti <-> ti, i.e., around the dash-dotted bisector. The function 
floor is the accidental triple rate _R 3 . 




FIG. 5: The expected experimental triple-coincidence rate 
N aa i(ti,t2,ti), obtained from the ideal triple-coincidence rate 
Pssi{t\,t2,ti) of Fig(4] by means of a moving bidimensional 
averaging window of area 4r,? and transition regions of size 
2rd- The central structure is in practice not related to the 
central peak of Fig(4] The scales of time axes in Fig(4] and [5] 
differ by more than three orders of magnitude. The floor of 
the function is again the accidental triple rate R 3 . 



N ssi (r) 
and N S i(r) 



R 3 (l + 2X) for short delays, 
R 3 (l + X) for long delays, 

X) for short delays, 



R 2 (l 
R 2 



for long delays, 



(8) 
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is therefore 
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4i? r c for short delays 
for long delays. 



(10) 



where "long" and "short" delays mean approximately 
|t| > t c + Td and \t\ < t c — respectively. It is 
not surprising that only the accidental-coincidence fac- 
tor, R 3 , and the signal-idler coincidence-peak factor, 
X = (2Rt c )^ 1 , show up in the formula for N SS i: all 
signal-signal effects are in fact too small to be detectable 
unless explicitly looked for. The graph of g c (t), EqJTl 



This conclusion is quantitatively confirmed by visual in- 
spection of [2, Fig. 3 and Fig. 4] (in the limit r c » t c l 
of coursefl]). Therefore, far from addressing the proba- 
bility of detecting two signal photons conditioned on the 
detection of an idler photons, these authors are again 
just measuring the temporal correlation between the sig- 
nal and idler arm. 
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Meaning of the conditioned second-order coherence 

Let us come back to the definition of second-order sig- 
nal coherence conditioned on the detection of an idler 
photon at time ij, given in Eq[3J This is the standard 
definition of Glauber's g^ 2 \ but over a post-measurement 
state obtained by applying the idler measurement oper- 
ator to the unconditioned SPDC density matrix: 



Pu 



tc\pEi(U)Ei(ti)] 



(11) 



The resulting state is of course not time invariant, be- 
cause of the introduction of the parameter tj. Equiva- 
lently, one can use modified observables for the evalua- 
tion of expectation values on the original state: 



(Y) P n 



(Ej(U)E t (U)) ' 



(12) 



The usefulness of the g^ 2 \ti,t 2 ) function for a time- 
invariant field is that it is insensitive to detector efficien- 
cies and depends only on r = t 2 — 1±, proportional to the 
probability (or rate, for a continuous field) of finding two 
photons at times t\ and t 2 . If <7^ (r) decreases (increases) 
in some neighbourhood of r = the field is said to be 
bunched (anti-bunched), meaning that photons "prefer" 
to come together (alone). If g^(0) is smaller (larger) 
than one the field statistics in a sufficiently small time 
window is sub-Poissonian (super-Poissonian) . 

However, due to the introduction of the idler time ti, 
the conditioned SPDC field is no more time invariant, 

(2) • r 

and g c turns out not to be proportional to the rate of 
arrival of signal photons at time t\ and t 2 , because the 
normalisation factors in the denominator of EqJ3] are not 
constant. For instance, if t 2 is set to a very large time, 
using EqsHJ[5J and [1] it is immediate to find that 



(13) 



irrespectively of the delay t\ — ti, although the detec- 
tion rate on the first signal detector is much larger when 



t\ = ti with respect to the case ti — ti 3> At, due to 

(2) 

signal-idler correlation. Therefore gc is not a good sub- 
stitute for detection rates; one should instead directly 
use the probability density of detecting signal photons at 
times t\ and t 2 provided an idler photon was detected at 
time ti, which, due to Bayes' theorem, is just the triple- 
coincidence rate P SS i(t\, t 2 ,ti) divided by R. 

It is also to be remarked that the last two special cases 
described in 0, pag.3, right], which should presumably 

(2) 

corroborate the idea that gc is a sensible figure of merit, 
can also be obtained from P SS i. Specifically, these exam- 
ples concern the fact that At is the signal-idler "coher- 
ence time" , 



P s si (- CO, (X)\U) 



R? 



and that unconditioned signal-signal events are thermal, 



Pssi{t\,t 2 \ - OO) 

P ss j(ti,oo| - oo) 



= 1 



R 2 {h -t 2 ) 
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= l + \gW( tl -t 2 )\ 



On the other hand, instead of showing any sign 
of sub-Poissonian behaviour, as the last special case 

f 21 

g c (ti,ti\ti) ~ seems to suggest, the point t\ = t% = t% 
is the maximum of P SS i- This is in agreement with the 
fact that one-arm SPDC radiation is thermal: if an idler 
photon was detected at time t{ and a signal photon at 
time t\ = U, then the most likely time for a second sig- 
nal photon to appear is just t 2 = ti. 



Conclusions 

It is shown that the experimental findings reported in 
2, Fig. 3 and 4] about conditioned signal-signal coinci- 
dences in an SPDC source are indeed related not to the 
conditioned second-order signal-signal coherence, but to 
the much simpler (to measure) and stronger signal-idler 
correlation. The author acknowledges fruitful discussions 
on the subject with H. Hiibel and I. M. Herbauts. 
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(see Fig[5]here). The authors of 3J claim that the peak 
at the centre of 0, Fig. 5] represents the contribution of 
multiple photon pairs, but this is clearly wrong, since it is 
just an artifact of the averaging method. In fact, the peak 
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